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Unlike  component  mode  substitution  methods,  no  approximating  assumptions 
are  made.  Thus,  natural  frequencies  and  mode  shapes  for  the  finite  element 
model  employed  are  the  same  with  or  without  the  substructuring  algorithm. 
This  is  demonstrated  by  computing  first  ten  natural  frequencies  and  the 
corresponding  mode  shapes  for  an  open  truss  helicopter  tail-boom  structure. 
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I. 


INTRODUCTION 


Analysis  of  structures  for  transient  dynamic  loads  is  of  considerable 
importance  in  many  fields  of  engineering.  Design  of  aircraft  structures, 
spacecraft  structures,  ship  structures,  machine  components  and  framed 
structures  for  earthquake  loading  are  some  examples  where  transient 
dynamic  analysis  procedures  are  required.  Many  of  these  structures  must 
be  represented  by  large  finite  element  models.  Therefore  it  is  desirable 
to  divide  the  structure  into  smaller  components  or  substructures  for 
efficient  dynamic  analysis. 

This  paper  addresses  the  question  of  dynamic  analysis  of  large 
structures  with  substructures  or  components.  A simple  and  efficient 
numerical  method  for  such  analysis  is  developed.  The  method  does  not 
require  generation  and  storage  of  stiffness  and  mass  matrices  for  the 
entire  structure.  It  requires  generation  of  only  substructural  stiffness 
matrices  and  storage  of  their  decomposed  factors.  Thus  the  method  is 
quite  compact  requiring  minimum  computer  storage. 

The  method  of  modal  superposition  is  adopted  for  dynamic  structural 
analysis.  In  this  approach  free  vibration  modes  and  the  corresponding 
natural  frequencies  must  be  first  computed.  This  information  is  then 
used  to  define  uncoupled  equations  of  motion  for  the  system.  Uncoupled 
equations  of  motion  are  then  easily  solved  using  Duhamel  Integral.  Thus 
the  problem  of  structural  dynamics  is  reduced  to  finding  natural 
frequencies  and  mode  shapes  with  substructures. 

Methods  of  component  mode  substitution  for  calculating  system 
eigenvalues  and  corresponding  eigenvectors  are  briefly  reviewed  in  the 
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next  subsection.  Section  II  defines  the  dynamic  structural  analysis 
problem.  the  method  of  dynamic  structural  analysis  with  substructures 
is  developed  in  Section  III.  An  example  problem  is  considered  in 
Section  IV  to  demonstrate  applicability  of  the  method. 

Component  Mode  Substitution  Methods 

In  component  mode  substitution  methods  of  free  vibration  analysis, 
the  structure  is  divided  into  a number  of  smaller  structures.  Mode 
shapes  of  each  substructure  are  calculated  for  a variety  of  boundary 
conditions.  Compatibility  and  equilibrium  conditions  at  the  inferfaces 
of  all  substructures  are  invoked  and  a reduced  eigenvalue  problem  for 
the  entire  structure  is  generated.  After  the  reduced  eigenvalue  problem 
has  been  solved,  system  mode  shapes  are  generated  using  previously 
ctefined  transformations  in  terms  of  component  modes. 

Advantage  of  component  mode  substitution  methods  is  in  the  experi- 
mental determination  of  system  mode  shapes.  It  is  easier  to  measure 
mode  shapes  and  natural  frequencies  of  smaller  components  as  compared  to 
measuring  mode  shapes  for  the  entire  system.  Knowing  mode  shapes  of 
each  component,  a method  is  needed  to  synthesize  the  system  modes.  The 
component  mode  substitution  methods  [1-7]  then  provide  such  a means  for 
generating  system  modes. 

Many  variations  of  the  component  mode  substitution  method  have  been 
presented  since  Hurty's  1960  paper  on  the  subject  [1].  Several  review 
articles  on  the  subject  have  appeared  in  the  literature  [2, 3, A, 5].  Some 
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other  recent  articles  that  deal  with  component  mode  substitution  concept 
are  the  Refs.  6 and  7. 

The  method  developed  in  the  paper  does  not  use  component  mode  substi- 
tution concept.  Therefore  details  of  these  techniques  are  not  presented 
here.  Reader  may  consult  Refs.  1-7  and  articles  cited  in  them  for  more 
details  on  component  mode  substitution  methods.  In  the  method  presented 
in  the  paper  one  calculates  eigenmodes  and  the  corresponding  eigenvalues 
for  the  system  with  the  substructural  data,  but  without  calculating  mode 
shapes  of  any  substructure.  No  assumptions  other  than  the  usual  assumptions 
of  linear  elastic  systems  modeled  by  finite  elements,  are  made. 
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II.  PROBLEM  OF  STRUCTURAL  DYNAMICS 

The  finite  element  technique  of  analysis  is  used.  The  equation  of 
motion  for  a typical  finite  element  model  for  the  structural  system  with 
n independent  degrees-of-f reedom  can  be  written  as 

Mu  + Cu  + Ku  = F(t)  (1) 

with  initial  conditions  specified  on  u and  u.  Here  M,  C and  K are 
(n  x n)  mass,  damping  and  stiffness  matrices,  respectively;  u,  u and  u 
are  the  displacement,  velocity  and  acceleration  vectors,  respectively 
and  F(t)  is  a forcing  function.  These  equations  are  derived  using  the 
Lagrangian  approach  [8],  If  it  is  assumed  that  the  system  described  by 
Eq.  (1)  is  linear,  then  M,  C and  K are  constant;  that  is  they  are  not 
functions  of  time  or  motion. 

Most  commonly  used  procedure  of  dynamic  analysis  of  linear  systems 
is  the  normal  modes  method.  This  is  also  referred  to  as  the  modal  or 
mode  superposition  analysis  and  is  adopted  in  the  present  paper.  In 
this  method,  one  first  considers  the  free  vibration  problem: 

Mu  + Ku  = 0 (2) 

Substituting 

u = 4>  sinw(t  - t ) (3) 

o 

into  Eq.  (2),  the  following  generalized  eigenvalue  problem  is  obtained: 

2 

K4>  * a)  M$ 


(4) 


The  n eigenvalues  of  Eq.  (4)  give  the  natural  frequencies  of  the  system 
and  the  corresponding  eigenvectors  are  the  mode  shapes.  The  complete 
solution  to  Eq.  (4)  can  be  written  as 
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2 

K<1>  = M4>ft  (5) 

i 2 2 

in  which  the  columns  in  4>  are  the  eigenvectors  4>  and  ft  = diag  (ok)  . 

The  eigenvectors  of  Eq.  (4)  are  linearly  independent  and  can  be 
orthonormalized  with  respect  to  the  mass  matrix  M.  These  vectors,  there- 
fore, form  a basis  for  the  n dimensional  vector  space.  The  vector  u(t) 
can  be  expressed  as  a linear  combination  of  the  eigenvectors  1 . This 
relationship  may  be  expressed  as  the  following  linear  transformation 

u(t)  = $n(t)  (6) 

where  n(t)eRn  is  usually  called  the  normal  coordinate  vector.  Substituting 

T 

Eq.  (6)  into  Eq.  (1)  and  pre-multipying  by  4>  , then  givs 

n + C*f|  + ft^n  = F*(t)  (7) 

T T 2 T 

Here  F*(t)  = 4>  F(t),  $ K4>  = ft  and  C*  = $ C$.  If  the  damping  matrix  C 
for  the  system  is  of  a restricted  form  [9,10],  then  C*  is  a diagonal 
matrix.  Equation  (7)  then  represents  a system  of  n uncoupled  equations. 
These  equations  can  be  solved  exactly  using  the  Duhamel  Integral  [8,10]. 

It  is  noted  that  the  finite  element  model  for  practical  systems  is 
generally  quite  large;  that  is  n is  generally  very  large.  For  many 
dynamic  analysis  problems,  it  is  reasonable  to  assume  that  only  a first 
few  eigenmodes  contribute  significantly  to  the  dynamic  response.  Thus, 
approximate  analyses  are  carried  out  with  only  p <<  n eigenvectors.  In 
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this  case,  4>  forms  an  M-orthonormal  basis  for  the  p dimensional  vector 
space  and  Eq.  (7)  represents  a system  of  only  p uncoupled  equations. 
Integration  of  Eq.  (7)  determines  the  normal  coordinates  peR^,  and 
Eq.  (6)  then  determines  the  real  displacement  u. 


III.  DYNAMIC  ANALYSIS  WITH  SUBSTRUCTURES 


The  problem  of  dynamic  structural  analysis  is  reduced  to  formulation 
of  the  eigenvalue  problem  of  Eq.  (4)  and  generation  of  the  desired  eigen- 
vectors and  eigenvalues.  Once  eigenvalues  and  eigenvectors  are  known, 
the  uncoupled  system  of  Eqs.  (7)  is  readily  obtained. 

Several  methods  for  solving  the  eigenvalue  problem  of  Eq.  (4) 
have  been  developed  in  the  literature  [11-15].  These  methods  take 
advantage  of  bandedness  property  of  the  K and  M matrices.  It  has  been 
recognized  that  it  is  desirable  to  divide  a large  structure  into  a 
number  of  smaller  structures  for  the  purpose  of  calculating  system 
eigenmodes  and  eigenvalues.  The  stiffness  and  mass  properties  of  a 
substructure  are  used  to  calculate  its  eigenmodes.  This  requires  solution 
of  several  smaller  eigenvalue  problems.  Eigenmodes  of  all  substructures 
are  then  synthesized  to  generate  system  eigenmodes.  Different  methods 
of  synthesis  generate  various  component  mode  substitution  methods  [5]. 

Most  of  these  methods  generate  approximate  eigenvalues  and  eigenvectors 
for  the  system. 

The  method  presented  in  the  paper,  however,  does  not  use  the  com- 
ponent mode  synthesis  concept,  since  the  eigenmodes  for  substructures 
are  never  computed.  The  method  directly  generates  the  desired  number  of 
system  eigenmodes  and  eigenvalues  using  the  substructural  stiffness 
matrices  and  the  mass  matrix  for  each  finite  element.  The  stiffness  or 
mass  matrix  for  the  entire  system  is  never  generated.  Moreover  the  mass 
matrix  for  any  substructure  is  not  required  to  be  generated  and  stored. 

The  method  developed  here  uses  the  subspace  iteration  algorithm  of 
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Ref.  12  for  solving  the  eigenvalue  problem  of  Eq.  (5).  The  basic  idea 
of  the  subspace  iteration  method  is  quite  similar  to  the  general  Rayleigh- 
Ritz  method.  In  the  method,  one  starts  with  a set  of  p « n linearly 
independent  vectors.  These  vectors  form  a basis  for  the  p-dimensional 
vector  space.  The  mode  shapes  for  the  system  are  expressed  as  linear 
combinations  of  the  assumed  basis  vectors.  The  parameters  of  these 
linear  combinations  are  obtained  by  solving  a reduced  eigenvalue  problem. 
An  iterative  process  is  then  used  to  improve  the  assumed  basis  vectors 
until  convergence  is  obtained. 

It  is  suggested  [12]  that  when  p mode  shapes  are  required,  one  should 
start  with  q > p linearly  indpendent  vectors,  where  q is  given  as 


q = min{2p,  p+8,  n} 


(8) 


This  improves  accuracy  of  the  first  p eigenvectors  and  the  corresponding 
eigenvalues.  Let  us  define  a linear  transformation 


X(0)$* 


(9) 


where  is  an  (n  x q)  matrix  whose  columns  are  the  assumed  basis 

vectors  for  R1"*  and  4>*  is  a (q  x q)  matrix  whose  elements  are  the  unknown 

multipliers  of  the  transformation.  Substituting  Eq.  (9)  into  Eq.  (5) 

(0)T 

and  premultiplying  by  X one  obtains 


K*4>*  = M*4>*fi* 


2 


(10) 


where 


K*  = X(0)  KX(0) 


M*  = X(0)  MX(0) 


(ID 


t 


The  superscript  indicates  the  iteration  number. 


j v5KT'  ~-7r 
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Equation  (10)  is  a reduced  eigenproblem  of  dimension  q that  is  solved  for  <t* 
2 

and  £1*  . Equation  (9)  is  then  used  to  obtain  approximate  system  mode 
shapes.  An  improved  set  of  basis  vectors  is  now  determined  from 


KX(1)  = Y(0) 


(12) 


where 


Y(°>  = M*<0) 


\ 

(13) 


The  iterative  process  is  continued  until  the  eigenvalues  determined 
from  Eq.  (10)  converge  to  within  a prescribed  tolerance. 


Subspace  Iteration  with  Substructuring 

In  using  the  subspace  iteration  algorithm  outlined  in  the  preceding 
with  substructuring,  one  needs  to  modify  calculations  of  Eqs.  (12)  and 
(13)  only.  If  one  can  perform  these  calculations  with  the  substructural 
matrices  rather  than  the  system  matrices,  then  one  has  a method  of 
structural  dynamics  with  substructuring.  Calculations  of  Eqs.  (12)  and 
(13)  can  indeed  be  performed  with  only  the  substructural  data.  This  is 
explained  in  the  following. 

Le  us  consider  Eq.  (13)  first.  Since  the  mass  matrix  M is  generated 
by  considering  contribution  to  it  from  each  finite  element  [10],  Eq.  (13) 
may  be  written  as 


Y 


(0) 


* 


(0) 


(14) 
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where  e Is  the  number  of  finite  elements,  m^  is  the  mass  matrix  for 
the  ith  element  and  A*  is  a Boolean  transformation  matrix  taking  contri- 
butions from  m*  into  M.  Note  that  for  each  i,  the  matrix  under  the 
summation  sign  is  an  (n  x n)  matrix.  Equation  (14)  may  be  rearranged 
as 


Thus  calculations  for  Y^^  from  Eq.  (15)  can  be  performed  without  ever 
generating  the  system  mass  matrix  M.  It  is  noted,  however,  that  calcu- 
lations of  Eq.  (13)  must  be  performed  in  each  subspace  iteration.  There- 
fore computational  effort  with  Eq.  (15)  is  greater  as  compared  to 
calculations  with  a system  mass  matrix  that  is  calculated  once  and 
stored  in  the  computer  core.  This  is  due  to  the  fact  that  calculations 
with  Eq.  (15)  essentially  imply  generation  of  M at  each  subspace 
iteration.  Thus,  there  is  a trade-off  between  the  computational  effort 
and  the  computer  core  requirement  for  calculations  of  Eq.  (13).  Note  that 
there  is  also  a possibility  of  computing  M only  once  and  storing  it  on 
some  external  device  such  as  a disk  or  magnetic  tape.  In  every  subspace 
iteration  the  stored  M can  be  readily  used.  This  procedure,  however, 
will  also  require  more  computational  time  as  compared  to  storage  of  M in  the 
core  because  of  the  time  delay  in  linking  and  rewinding  the  external  device. 

Calculations  for  from  Eq.  (12)  can  also  be  performed  with  only 

substructural  matrices.  The  reason  for  this  is  that  one  can  view 
Eq.  (12)  as  an  equilibrium  equation  governing  static  response  of  the 
system  to  a pseudo-load  matrix  Y^^  whose  each  column  represents  a 
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loading  conditions.  Equilibrium  equations  of  the  type  of  Eq.  (12)  can 
be  solved  using  the  substructuring  concept  [10].  The  reader  who  is 
familiar  with  static  structural  analysis  with  substructures  will  readily 
recognize  the  following  development.  Others  may  consult  Ref.  10  for 
more  discussion  on  static  structural  analysis  with  substructures. 

For  the  purpose  of  solving  from  Eq.  (12)  using  the  substructuring 

concept,  one  divides  the  entire  structure  into  several  substructures. 

Each  substructure  is  viewed  as  a hyper-finite  element  that  is  connected 
to  other  substructures  only  at  its  interfaces.  For  each  substructure  a 
set  of  interior  generalized  coordinates  and  another  set  of  interface  or 
boundary  generalized  coordinates  is  defined.  The  hyper-element  force- 
displacement  relationships  are  derived  in  terms  of  only  the  boundary 
generalized  coordinates  by  using  essentially  the  static  condensation 
procedure  [10].  The  hyper-element  force-displacement  relationships  are 
synthesized  to  generate  a reduced  system  of  equilibrium  equations  in  terms 
of  only  the  boundary  generalized  coordinates  for  the  system.  This  reduced 
system  of  equations  is  solved  for  boundary  generalized  coordinates  and 
following  the  usual  finite  element  approach,  displacements  at  all  points 
of  the  structure  are  calculated.  This  procedure  is  developed  in  the 
following. 

First,  let  us  partition  Eq.  (12)  as  follows: 


(16) 


4m 


* 
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where  subscripts  B and  1 refer  to  variables  associated  with  the  boundary 
and  interior  generalized  coordinates  for  the  entire  structure.  Now, 
using  the  second  line  of  Eq.  (16),  is  eliminated  from  the  first 

line  to  obtain 


v?5  - fe 


(17) 


where 


fb = ybo) + <*mo) 


kb  - kbb  + SxQ 


- K_1  K 
II  kib 


(18) 


(19) 


(20) 


It  is  interesting  to  note  that  matrices  Fg  and  in  Eqs.  (18)  and  (19) 
can  be  generated  using  the  substructural  matrices  without  ever  generating 
FBB’  Kn  or  ^BI  For  entlre  structure.  For  this  purpose  one  also 
partitions  the  equilibrium  equation  for  the  rth  substructures  as 


^B 

■4i' 

rY(Dr1 

ab 

rY(°)r- 

B 

5b 

KII_ 

y(l)r 

AI 

Y(°)r 

_I 

(21) 


where  superscript  r refers  to  the  rth  substructure.  From  the  second  line 
of  Eq.  (21),  one  obtains 
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y 


Y(1>r  _ Y(°)r  VT  y(Dr 

K1IXI  ■ Yi  ' kib  x 


(22) 


(1)  r 

Eliminating  from  the  first  line  of  Eq.  (21),  using  Eq.  (22),  one 

obtains 


KB  *B  _ FB 


(23) 


where 


KB  = ^B  + 


(24) 


FB  = YB°)r  + Yl°)r 


(25) 


* ' "\ 


Qr  " - lKu1_1  4 


(26) 


Equation  (23)  is  a force-displacement  relationship  for  the  rth  sub- 
structure in  terms  of  its  boundary  generalized  coordinate.  The  boundary 
stiffness  matrix  and  the  effective  boundary  load  matrix  Fg  for  the 
entire  structure  are  synthesized  from  the  substructural  data  of  Eqs. 

(24)  and  (25)  using  the  standard  finite  element  procedures: 


L T 

\ sr 

r=l 


(27) 


and 


F»  ■ ,b0)  + X “rT  Q'  v50)r 


< 28) 


:-r  • 


sr  -•  yawn  < 
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Is 


where  L is  the  number  of  substructures  and  6 is  a Boolean  transformation 

r 

matrix.  It  is  noted  here  that  the  matrix  in  Eq.  (26)  is  usually  not 

inverted.  This  matrix  is  symmetric  and  banded.  Efficient  procedures 

r r 

are  used  to  decompose  and  to  solve  for  Q in  Eq.  (26)  using  forward 

and  backward  substitution.  Also  is  usually  quite  sparse  and 

IB 

advantage  of  this  property  is  realized  in  calculations. 


Now 


r ( 1 ) 


the  matrix  Xg  is  calculated  from  Eq.  (17).  The  matrices 


(1  \ r 

Xj  , r = 1,  2,  ...»  L,  are  calculated  from  Eq.  (22).  Finally  the 

matrix  is  assembled  from  X^^  and  X^^r.  It  is  noted  that  the 

r r r r 

matrices,  such  as  decompositions  of  K and  K , Q , It,  and  6 are 

B II  IB 

available  from  the  static  structural  analysis  for  use  in  the  preceding 
equations.  Thus,  the  proposed  subspace  iteration  algorithm  with  sub- 
structuring blends  quite  nicely  into  the  static  structural  analysis 
with  substructuring. 

The  subspace  iteration  algorithm  with  substructures  is  then 
summarized  as  follows: 


Step  1 Start  with  an  (n  x q)  matrix  X^  whose  columns  are  estimates 
of  q eigenvectors. 

Step  2 Compute 


*<°>  - M*<°> 


(29) 


using  finite  element  operations  similar  to  that  of  Eq.  (15). 


Step  3 Solve  for  from  the  equation 


IS 


T’ 


i E 


y ' 


. ‘i-v' 


5“  * - 
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KX(1)  - Y(0> 


(30) 


using  Eqs.  (17)  and  (22), 


Step  4 Compute 


Y(1)  = MX(1) 


(31) 


using  finite  element  operations  of  Eq.  (15), 


Step  5 Calculate  the  following  (q  x q)  matrices 


K*-XW>  Y(0) 


M*  = X(1^  Y^ 


(32) 


Step  6 Solve  for  all  eigenvalues  and  eigenvectors  of  the  projected 
eigenvalue  problem  of  Eq.  (10).  Note  that  the  generalized 
Jacobi  iteration  [14]  or  any  other  method  of  solving  an  eigen- 
value problem  may  be  used  in  this  step. 

Step  7 Compute 


x(!)  = ^(1)  ^ y(l)  = y(l) 


(33) 


Step  8 Check  for  convergence  of  eigenvalues.  If  all  the  ratios 


,(D 


,(0) 


- in, 


,(D 


(34) 


are  within  a specified  tolerance,  then  stop  the  iterative 
process.  Otherwise  return  to  Step  3 with  Y^  = Y^. 
After  convergence,  the  first  p columns  of  are  the 


:U 


- — I,. — 


— 


'**v' 


_____ 
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required  eigenvectors  and  the  first  p elements  in  ft*  are  the 

f 

corresponding  eigenvalues. 

The  rate  of  convergence  of  the  subspace  iteration  algorithm  depends 
on  how  close  the  starting  subspace  is  to  the  eigenvectors  of  the  system. 
Many  times  eigenvectors  from  previous  analyses  are  available  and  form  an 
excellent  starting  subspace.  For  many  structural  problems,  authors  have 
simply  used  first  q columns  of  an  (n  x n)  identity  matrix  as  the  starting 
subspace  [16-18].  The  procedure  suggested  in  Ref.  12  had  also  been  used 
successfully  [19,20],  Here  the  first  column  of  contains  simply 

diagonal  elements  of  the  mass  matrix  and  other  columns  are  unit  vectors 
with  +1  at  the  coordinates  with  largest  ratios  ra^/k^. 
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IV.  AN  EXAMPLE  PROBLEM 

As  an  example  of  calculating  eigenvalues  and  eigenvectors  with  sub- 
structures, an  open  helicopter  tail-boom  structure  is  considered.  The 
structure  is  modeled  by  108  truss  members  and  28  joints  as  shown  in 
Figure  1.  The  structure  has  72  degrees  of  freedom,  so  the  eigenvalue 

problem  of  Eq.  (4)  is  of  dimension  72  * 72.  For  sample  calculations, 

2 

cross-sectional  area  of  all  members  is  1.0  in.  , the  material  weight 

3 6 

density  is  0.10  lb/in.  and  Young's  modulus  is  10  psi. 

The  first  10  natural  frequencies  and  the  corresponding  mode  shapes 
for  the  structure  are  obtained  using  the  subspace  iteration  algorithms 
with  and  without  substructuring.  Two  computer  programs  based  on  the 
subspace  iteration  have  been  developed,  one  without  substructui ing  and 
the  other  with  substructuring.  These  programs  are  used  to  obtain  the 
results  reported  herein.  Note  that  the  generalized  Jacobi  iteration  [14] 
is  used  to  solve  the  eigenproblem  of  Eq.  (10)  in  both  programs. 

For  the  algorithm  without  substructuring,  the  half-bandwidth  of 
matrices  K and  M is  21.  Thus  72  * 21  stiffness  matrix  is  computed.  It 
is  decomposed  and  stored  for  use  in  Step  3 of  the  algorithm.  The  72  * 21 
mass  matrix  for  the  structure  is  also  computed  and  stored  for  use  in 
Steps  2 and  4 of  the  algorithm. 

For  the  algorithm  with  substructuring,  the  tail-boom  structure  is 
divided  into  three  substructures  by  partitioning  it  at  nodes  9-12  and 
17-20  as  shown  in  Fig.  2.  Nodes  25-28  are  also  treated  as  boundary  nodes. 
Figure  2 shows  global  and  local  (or  substructural)  numbering  systems  for 
various  substructures.  Substructure  1 has  36  members,  12  boundary  degrees 
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Finite  Element  Model  for  the  Tail-Boom  Structure 
(Degrees  of  freedom  for  joints  5 and  28  are 
1*  2,  3 and  70,  71,  72,  respectively.) 


Figure  1.  Helicopter  Tail-Boom  Structure 


(a)  AN  OVERALL  NUMBERING  SYSTEM  FOR 
BOUNDARY  NODES. 


SUBSTRUCTURE  1 


SUBSTRUCTURE  2 SUBSTRUCTURE  3 


(b)  NUMBERING  SYSTEM  FOR  BOUNDARY  AND 
INTERIOR  NODES  FOR  EACH  SUBSTRUCTURE. 


NOTE:  FOR  CLARITY  DIAGONAL  MEMBERS  ARE  NOT  SHOWN. 


Figure  2.  Nodal  Numbering  Systems  for  Substructural  Formulation 
for  the  Finite  Element  Model  of  the  Helicopter 
Tail-Boom  Structure 


20 


of  freedom  and  24  interior  degrees  of  freedom.  Substructures  2 and  3 

each  have  24  boundary  degrees  of  freedom,  12  interior  degrees  of  freedom 

and  36  members.  The  boundary  stiffness  matrix  K for  the  entire  struc- 

B 

ture  is  a 36  x 36  matrix  with  half-bandwidth  of  12.  Decomposed  part  of 
this  matrix  is  stored  for  use  in  Eq.  (17).  The  decomposed  part  of  the 
matrix  for  each  substructure  is  also  stored  for  use  in  Eq.  (26). 

r 

The  matrix  Q is  calculated  from  Eq.  (26)  and  stored  for  each  substruc- 
ture. In  this  program  the  mass  matrix  for  the  structure  is  not  stored. 
Rather,  calculations  of  Eqs.  (29)  and  (31)  are  carried  out  memberwise. 
Therefore,  one  should  expect  computational  time  to  be  slightly  greater 
with  this  program. 

The  starting  subspace  for  both  programs  is  obtained  by  taking  the 
first  18  columns  of  a 72  * 72  identity  matrix.  It  should  be  noted,  how- 
ever, that  degrees  of  freedom  with  substructuring  formulation  are  numbered 
differently  as  compared  to  formulation  without  substructuring.  Therefore, 
the  starting  subspace  for  the  two  programs  is  different.  However,  the 
two  programs  converged  to  exactly  the  same  natural  frequencies  and  the 
mode  shapes.  The  results  are  summarized  in  the  Table.  The  ten  eigen- 
vectors are  listed  in  the  Appendix.  The  tolerance  limit  for  eigenvalue 
convergence  in  Step  8 was  10  ^ for  both  programs. 

For  the  results  given  in  the  table  the  subspace  iteration  algorithm 
with  substructuring  took  9 iterations  to  converge.  The  computing  time 
was  72.3  sec  on  IBM  3 7 0-1 68(G)  computer.  The  program  based  on  the  algo- 
rithm without  substructuring  took  13  iterations  to  converge  with  a com- 
putational time  of  82.5  sec.  The  reason  for  four  more  iterations  with 
this  program  is  that  the  starting  subspace  is  apparently  not  as  good  as 


TABLE  NATURAL  FREQUENCIES  WITH  OR  WITHOU' 


SUBSTRUCTURE 


Natural  Natural 

Frequency  Frequency  in  Hz 

No . 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


21.8236 

23.1330 

101.5989 

105.2161 

107.3932 

200.8760 

227.5787 

239.3666 

241.9844 

377.4557 


----- 





for  the  previous  program.  However  when  the  same  starting  subspace  is 
chosen  for  the  second  program,  the  algorithm  converges  to  the  same  solu- 
tion in  9 iterations  with  a computational  time  of  61.0  sec.  This  time 
is  less  than  the  computing  time  with  the  substructuring  program.  The 
reason  is  that  with  the  substructuring  program,  the  mass  matrix  is  cal- 
culated in  each  subspace  iteration  for  use  in  Step  4 of  the  algorithm. 
When  the  substructuring  program  is  modified  to  calculation  of  the  mass 
matrix  only  once,  the  computing  time  is  reduced  to  57.8  sec. 
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V.  DISCUSSIONS  AND  CONCLUSIONS 

An  efficient  numerical  procedure  for  structural  dynamics  with  sub- 
structuring is  developed.  The  method  is  based  on  the  subspace  iteration 
algorithm  for  calculating  eigenvalues  and  eigenvectors  of  a general 
eigenproblem  of  Eq.  (4).  No  approximating  assumptions  are  made  in 
developing  the  procedure.  Thus,  for  a given  finite  element  model  of 
the  structure,  the  natural  frequencies  and  the  corresponding  mode  shapes 
are  the  same  with  or  without  the  substructuring  method. 

The  method  proposed  herein  does  not  use  the  component  mode  sub- 
stitution ideas  developed  during  the  past  18  years.  Thus,  solution  of 
the  eigenproblem  for  each  component  is  not  necessary. 

Conceptually  the  method  is  fairly  straightforward  and  simple  to 
program.  Another  advantage  of  the  method  is  that  it  uses  most  of  the 
data  that  is  already  available  from  structural  analysis  for  static 
response.  Thus  the  method  blends  quite  nicely  into  static  structural 
analysis  with  substructures.  On  the  other  hand  this  will  not  be  true 
if  component  mode  substitution  methods  are  used,  as  some  additional 
calculations  with  substructures  would  be  necessary  for  obtaining  static 
response  of  the  structure. 

It  should  be  noted  that  the  method  of  calculating  natural  frequen- 
cies and  eigenvectors  with  substructures  can  be  easily  integrated  into 
the  optimal  design  algorithm  with  substructures  [21]  to  treat  constraints 
on  natural  frequencies  of  the  structure.  Also  optimal  design  procedures 
with  substructuring  can  be  now  developed  for  structures  subjected  to 
• transient  dynamic  loads. 
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